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Abstract. 

The dynamical state of the solar nebula depends critically upon whether or not 
the gas is magnetically coupled. The presence of a subthermal field will cause laminar 
flow to break down into turbulence. Magnetic coupling, in turn, depends upon the 
►^. . ionization fraction of the gas. The inner most region of the nebula (< 0.1 AU) is 

r — ■ magnetically well-coupled, as is the outermost region (> 10 AU). The magnetic 

status of intermediate scales (~ 1 AU) is less certain, ft is plausible that there is a 
(f) ■ zone adjacent to the inner disk in which turbulent heating self-consistently maintains 

^»0 ' the requisite ionization levels. But the region adjacent to the active outer disk is 

<^^> likely to be magnetically "dead." Hall currents play a significant role in nebular 

0"N ■ magnetohydrodynamics. 

^^ ' Though still occasionally argued in the literature, there is simply no evidence to 

C| support the once standard claim that differential rotation in a Keplerian disk is prone 

O ,■ to break down into shear turbulence by nonlinear instabilities. There is abundant 

evidence — numerical, experimental, and analytic — in support of the stabilizing role 
of Coriolis forces. Hydrodynamical turbulence is almost certainly not a source of 
enhanced turbulence in the solar nebula, or in any other astrophysical accretion 
disk. 
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1. Introduction 

Despite decades of activity, a reliable model for the dynamical state 
of the primitive solar nebula remains elusive. It is of course essential 
to have such an understanding if we are to make sense of nebular 
chemistry transport processes, planet formation and migration, and 
evolutionary history. The nebular temperature profile is of particu- 
lar importance, since much of the meteoritic data (e.g. condensation 
properties of refractory elements) bear directly upon this. 

A significant source of heat for the solar nebula will be turbulent 
dissipation, at least in those regions of the disk prone to instability. 
This turbulence also results in angular momentum and outward energy 
transport at rates far in excess of those associated with molecular vis- 
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cosity; one speaks therefore of an enhanced turbulent, or "anomalous" 
viscosity. These properties are the essence of an accretion disk. They 
allow the central star and disk to evolve, and in some circumstances 
(e.g. FU Orionis outbursts) the evolution may be eruptive. In this con- 
tribution, we review developments in our understanding of accretion 
disk turbulence as it applies to the solar nebula, and by extension, 
to protostellar and protoplanetary disks more generally. Recent gen- 
eral disk reviews may be found in Papaloizou & Lin (1995), Lin & 
Papaloizou (1996), and Balbus & Hawley (1998). 



2. Accretion Disk Turbulence 

2.1. A Preliminary Comment 

For many years, it has been the custom for disk articles to identify the 
source of anomalous viscosity as the principal difficulty of accretion 
disk theory, use words like "mysterious" to describe its origin, and pay 
homage to nonlinear instabilities, convection, and possibly magnetic 
fields. The author would then move on to a simple viscous model, which 
would form the basis of the discussion. 

Substantial progress has been made in the last decade. The ori- 
gin of accretion disk turbulence is not mysterious, it is in fact very 
simple. The presence of any kind of weak magnetic field (in a sense 
described below) in a conducting, differentially rotating fluid renders 
it violently unstable. Remarkably, the instability was known at some 
level (though its consequences were certainly not appreciated) forty 
years ago (Velikhov 1959, Chandrasekhar 1960). Its importance to ac- 
cretion disks was noted by Balbus & Hawley (1991, 1992a,b), and its 
subsequent elucidation and development have been carried through by 
many investigators (e.g., Goodman and Xu 1994; Hawley, Gammie, and 
Balbus 1995; Brandenburg et al. 1995; Matsumoto et al. 1996; Stone 
et al. 1996; Gammie 1996; Terquem and Papaloizou 1996; Ogilvie and 
Pringle 1996; Balbus and Hawley 1998 for review). 

If earlier work has inspired confidence in magnetohydrodynamic 
(MHD) turbulence, then the most recent findings compel outright chau- 
vinism. Not only is it established that weakly magnetized disks become 
unstable, developing enhanced turbulent transport, it also clear that 
hydrodynamical turbulence is fundamentally unsuitable to serve in 
this role. One may as well try to excite convection in a top-heated 
fluid as excite shear turbulence in a Keplerian disk 1 . Neither ven- 



1 This pithy observation was apparently noted by Sir Harold Jeffreys in conver- 
sation with L. Mestel (Papaloizou, private communication). 



ture will work, "nonlinear instabilities" notwithstanding. This rather 
unequivocal position will be justified later in this paper. 

2.2. The Magnetorotational Instability 

We begin with a brief presentation of the linear magnetic ("magne- 
torotational") instability, almost certainly the basis for any model of 
turbulent transport in the nebula. Set up a standard cylindrical coor- 
dinate system (R, <p, z) centered upon a central gravitating body. The 
equations of motion for an element of fluid orbiting in a plane in the 
field of central potential <& are 

d 2 R /d<p\ 2 d$ „ 

d 2 ct> dRdcj) _ 

R ^ + 2 Htt- f ^ (2) 

where fn and fs represent local forces in the R and (j) directions. We 
allow the particle to make small excursions from a circular orbit R = 
i?0j = ^ot by introducing small quasicartesian x and y variables 

R = R + x 4) = n t + y/R . (3) 

Substituting equations (3) into equations (1) and (2), canceling lead- 
ing order terms, and retaining those linear in x and y leads to the Hill 
equations. These were first used to study the Sun-Earth-Moon triplet 
system: 

dn 2 

» + 2S!± = / r (5) 

We follow the convention of writing dots for time derivatives and have 
resolved the local force into its x and y components. The angular ve- 
locity gradient dfl 2 /dlnR refers to the circular velocity profile Q(R), 
and we have suppressed the "0" subscript for clarity. 

In what follows, fluid pressure forces are unimportant, and we need 
consider only magnetic stresses in the /-forces. Consider an equilibrium 
vertical magnetic field. If we restrict our attention to fluid displace- 
ments in the plane of the disk, varying only with z and t as e H wi - fcz ) ; 
these stresses take a very simple form: the magnetic tension force is 
— (k • ua) 2 $ (e.g., Balbus and Hawley 1992a), where ua is the Alfven 
velocity 

UA = ^r P > (6) 



and £ is the two-dimensional displacement vector. The equations of 
motion take the form 



x-2Vty 



dCt 2 



dlnR 
y + 2ttx = -(k • u A fy 



k • ua) 2 x, 



(7) 

(8) 



These are equations for classical Alfven waves (more precisely an Alfven 
and a slow wave, which is degenerate with the former), modified by 
Coriolis forces and a radial tidal field. The magnetic tension terms 
proportional to (k • ua) 2 act like spring forces, with a tunable spring 
constant depending upon the magnitude of k. 

Equations (7) and (8) may be satisfied provided that u satisfies 



U) 



u 2 [n 2 + (k-u A ) 2 ] + (k- u A ) 



(k ■ u A ) 2 + 



dn 2 



dlnR 



a simple quadratic in uj 2 . Here, n 2 is the epicyclic frequency, 



i dR A n 2 



K 



R 3 dR 



If 



(k • u A ) 2 + 



dn 2 

dlnR 



<0, 



0, (9) 



(10) 
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then the system will be unstable. Physically, the left side equation (11) 
is (with a minus sign) the net centrifugal force on a displaced element, 
the magnetic field here appearing as a restoring, stabilizing agent. The 
inequality is a statement that this force should have a radial outward 
component for destabilization. 

Assuming that large enough wavelengths are allowable (a weak spring) 
when the angular velocity (not the angular momentum) decreases out- 
wards, the disk will be unstable. But of course that is the behavior of a 
Keplerian (or any other astrophysical disk) rotation law. Only slightly 
more work is needed to reveal that there is a maximum growth rate 



1 dSl 



uj. 



2dlnR 
which, for a Keplerian rotation law Q ~ r ~ 3 / 2 ^ [ s obtained when 

A 2 



(k • u A ) 



16 



(12) 



(13) 



The growth rate (12) is quite large, amounting to an amplification 
factor of about 100 per orbit. 



In what sense must the field be "weak" for the instability to run? 
Equation (11) provides the key. The field must be sufficiently weak 
that the Alfven tension frequency is less than y/3fl in a Keplerian disk. 
Wavenumbers in excess of this limitation are stabilized by the tension 
forces. If the field is too strong, all wavelengths able to fit inside the thin 
disk will become stabilized. Generally, this restricts the unperturbed 
Alfven speed to be less than of order the disk sound speed (Balbus and 
Hawley 1991). 

To understand the instability more fully, go back to our connecting 
spring analogy. When the spring is stretched so that one mass is orbiting 
slightly farther out than the other, the inner mass rotates more rapidly 
on the lower orbit. The tug of the spring pulls back on this mass, 
causing it to lose angular momentum. (This is the torque on the right 
hand side of equation [8]). This loss compels the mass to descend to 
a yet lower orbit, while the gain of angular momentum in the outer 
mass sends it to a higher orbit. The spring stretches more, and the 
process runs away. Note that this is an instability only if the spring is 
relatively weak; if the spring is too strong the mass points will simply 
oscillate as they orbit. The instability is strongest if the spring constant 
is comparable to f2 2 , as indicated by equation (13). Note as well that 
angular momentum transport is not some nonlinear outcome of the 
instability, it is the essence of the instability. 

Second, the instability is even more robust than our simplified pre- 
sentation indicates (Balbus and Hawley 1998). The maximum growth 
rate (12) is independent, not only of the strength of the magnetic field, 
but also of the geometry of the magnetic field. Even a purely toroidal 
field is unstable (Balbus and Hawley 1992b; Terquem & Papaloizou 
1996), with the same maximum growth rate (Balbus & Hawley 1998). 
The presence of buoyancy or other pressure effects does not affect the 
stability criterion, since the most unstable displacements lie in the 
plane of the disk. Furthermore, it has been conjectured that equation 
(12) represents the maximum possible growth rate that any instability 
tapping into the free energy of differential rotation can achieve (Balbus 
& Hawley 1992a). 

Finally, three-dimensional numerical simulations of the nonlinear 
phase of the instability (Hawley et al. 1995), combined with analytic 
arguments (Balbus and Papaloizou 1999) eliminate any doubt that the 
instability leads to turbulence, that the turbulence is fully developed, 
and that significant angular momentum transport results. Values of the 
"a parameter" (cf. next section) obtained in local simulations extend 
over the range 5 x 10~ 3 < a < 5 x 10 _1 (Brandenburg et al. 1995; 
Hawley et al. 1995), with lower values corresponding to initially uniform 
toroidal fields or fields with vanishing mean values, and the higher 
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values corresponding to the presence of a mean axial field. If the weak 
field MHD turbulence is present, it will act like a greatly enhanced 
disk viscosity, both diffusively and dissipatively. It will dictate the 
local thermal structure of the solar nebula. The question of importance 
becomes, where in the disk will the ionization levels be sufficiently high 
to permit good magnetic coupling? 

2.3. Turbulent Heating in the Nebula 

Viscous, turbulence-enhanced accretion has been a popular and stimu- 
lating model for the post-infall global evolution of the solar nebula (Lin 
and Papaloizou 1980). Classical viscous accretion disk models allow one 
to determine macroscopic disk properties in terms of two constants: the 
accretion rate M 2 , and a dimensionless disk viscosity parameter, a. 
The latter is related to the (turbulent) viscosity vt by 

vt = acsH (14) 

where cs is the isothermal sound speed, and H is the disk scale height 
(Pringle 1981). (The quantity cs is computed with a mean mass per 
particle appropriate to a molecular hydrogen gas of cosmic abundances, 
~ 2.33m p , where m p is a proton mass.) If, as we believe, MHD turbu- 
lence is the physics behind a, then it should be possible to go beyond 
simple parameterizations. 

We know, for example, that the disk must be hot enough to maintain 
an ion population capable of coupling the field to the molecular nebular 
gas. Where is this requirement self-consistently met? Using an a model 
as our starting point, let us parameterize our results in terms of the 
column density £ and midplane optical depth r. The midplane nebula 
temperature T may be found from energy conservation: 

4 SirR^a 

and angular momentum conservation, 

M = SnZacsH. (16) 

(Shakura and Sunyaev 1973; Frank, King, and Raine 1992). Here G is 
the gravitational constant, M is the central mass, and a is the Stefan- 
Boltzmann constant. Equation (15) expresses the condition that surface 
radiation losses are balanced by internal mechanical energy dissipation; 



2 Note that this accretion rate refers to radial drift through the disk, not an 
external infall rate. 
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equation (16) employs a vanishing stress boundary condition at the 
inner edge of the disk. 

These two equations lead to a temperature of 

T = n « K (^)" i (l)" ! * < 17) 

where Rau is the disk radial location in astronomical units. From 
the point of view of magnetic coupling, this is an interesting tem- 
perature regime, neither prohibitively cool, nor dominated by ions. 
Clearly, different regions of the nebula will have very different mag- 
netic properties, very different internal temperatures, and very different 
dynamical states. The solar nebula is likely to have been highly inho- 
mogeneous. Among other consequences, the ease with which a planet 
migrates through the nebula depends very much upon whether the gas 
is turbulent or not (Terquem, Papaloizou, and Nelson, this volume). 
With the database of new planets and even new planetary systems 
rapidly swelling, understanding the extent over which protoplanetary 
disks maintain turbulence is a key problem. 



3. Resistivity and Hall Currents 

Outside of the first few 0.1 AU, the ionization fraction in the solar 
nebula is likely to be very small. Of course very little ionization is 
needed to ensure good coupling in astrophysical plasmas, but even this 
may not always be attainable. There are two important effects in a low 
ionization plasma which bear discussion. 

3.1. Resistivity 

The resistivity of partially ionized plasma is (e.g., Blaes & Balbus 1994) 

r) = 230 (^] T 1 ' 2 cms" 2 , (18) 

where n n is the neutral number density, n e is the electron number 
density, and T is the temperature. A measure of the relative importance 
of the resistivity is given by the magnetic Reynolds number 

Hr r 2 

ReM = f^l^ £_, (19) 

where the final expression is our operational definition. When Reu is 
smaller than Cg/u 2 ^, resistivity affects the most rapidly growing linear 
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wavelengths, and when Reu falls below 47r 2 , the linear instability is 
completely quenched (Stone et al. 1999). (Note that c 2 s /u\ is generally, 
though not inevitably, larger than 47r 2 .) On scales of 1 AU in the solar 
nebula, this corresponds to ionization fraction of about 10~ 13 . [This 
is a small number, but it is large enough to be in a regime where we 
expect dust grains to be playing a relatively minor role (Umebayashi 
&; Nakano 1988)]. The dominant ions will be Na + and K + , abundant 
alkalis with low ionization potentials. 

The ionization in the inner nebula will be thermal [tenuous outer 
layers may be regulated by cosmic or X-rays, (Gammie 1996, Glassgold 
et al. 1997)], and an ionization fraction of 10~ 13 is associated with 
temperatures of 900-1100 K over a broad range of densities. Taking this 
result together with equation (17), we find that formally, an accretion 
disk with a = 10~ 2 can maintain this minimal ionization out to radii 
of order 1 AU. 

We need to be careful. Conditions for linear instability to be present 
and conditions for maintaining turbulence need not be the same. (One 
example of this is the linear two-dimensional magnetorotational in- 
stability, which cannot sustain turbulence because of the anti-dynamo 
theorem. Quite generally, shear turbulence requires three dimensions.) 
A local study by Fleming, Stone, & Hawley (1999), has shown a marked 
decrease in the level of MHD turbulence for R&m < 5 x 10 4 , with tur- 
bulence decaying once Reu < 2 x 10 4 . These values are far too high to 
affect the linear behavior of the instability, and indeed the simulations 
are indistinguishable early in the linear regime. It is possible that the 
upper limit represents the discreteness of the grid, since the behavior of 
R&M = 5x 10 4 and Li Reu = °o" (i.e., no explicit resistivity) are likewise 
indistinguishable. Thus, finite magnetic Reynolds number effects may 
be present at higher values as well. 

Reu is extremely sensitive to the temperature. If the ionization 
fraction is much less than the abundance of K (~ 10~ 7 ), the Saha 
equation gives 



,., |/o ; ! /'2.4 x 10 15x xrl 
6.47 x IO^VCTo 7 
n„ \ n 



exp(-25, 188/T) 



1.15 x 10 



-ii 



(20) 



where a_7 is the K abundance in unit of 10~ 7 , T3 is the temperature 
in units of 10 3 K; the Boltzmann factor is normalized to its value at 
1000 K, and the neutral density normalization is given by fundamental 
constants of nature which happen to lead to a value within an order of 
magnitude (or two) of that expected for the solar nebula. This leads to 



a magnetic Reynolds number of 

"" "exp(-25, 188/T) 



Re M = 15.8 T" 7 a_!:R 



5/4 1/2-3/2/ 2.4 xl0 15 \ 1/: ' 



3 U,_ 7 -LL AU ( 

'fcr. 



■ 



1.15 x 10 



-ii 



(21) 



At face value this is well below what the Fleming et al. simulations 
suggest is needed for unencumbered turbulence, but the exponential 
Boltzmann acts like a switch. One need not venture more than a few 
tenths of an AU inward for conditions to change dramatically, say, at 
T = 1400 K. We are then in the regime of full MHD turbulence. 

This behavior is the most salient difference between a simple a model 
for the nebula, and one based upon MHD turbulence: a very sharp 
change in the dynamical state of the inner 1 AU. Furthermore, it is also 
self-consistent to have a quiescent disk as we move inward of 1 AU, with 
much lower temperatures, and much lower ionization fractions. This 
would mean a smaller value for the turbulent stress, which is the basis of 
the self-consistency. More turbulence means more ions to help maintain 
the magnetic coupling, less turbulence means the converse. Clearly, 
the transition between quiescence and turbulence may be explosive: a 
slight increase in temperature turns into an increase in ionization, an 
increase in the alpha parameter, and yet higher temperatures. Gammie 
and Menou (1998) have argued that dwarf novae outbursts are related 
to just such a transition to MHD turbulence; it is very possible that 
FU Orionis outbursts are a similar phenomenon. 

3.2. The Hall Effect 

The second effect may be understood as follows. In the absence of 
resistivity, the standard induction equation for the magnetic field B is 

— = Vx(vxB), (22) 

where v is the fluid velocity. The question is, which fluid? Under most 
circumstances, the question is moot because the difference between 
ion, electron, and neutral velocities is too small to matter. But when 
the ionization fraction becomes sufficiently small (and 10~ 13 qualifies), 
we must be careful with the distinction. A more accurate induction 
equation is 

— = Vx(v e X B), (23) 

that is, the electron velocity replaces the "fluid" velocity. This is be- 
cause the magnetic field lines are tied most effectively to the most 
mobile charge carriers. In the case of interest, the inertia of the charge 
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carriers is negligible, and the fluid velocity is to a very good approxi- 
mation that of the neutrals. Since 

v e = v + (v e - Vi) + (vi -v), (24) 

the electron velocity may be represented as the sum of the neutral 
velocity, a second term proportional to the current {vi is the ion veloc- 
ity), plus a final term determined by the strength of ambipolar diffusion 
coupling. In interstellar molecular gas, this final term generally turns 
out to be more important than the second. But in protostellar disks, 
the second term, which gives rise to the Hall effect in the laboratory, 
is dominant. This point is made, using the formalism of conductivity 
tensors, in an important recent paper by Wardle (1999). We will shall 
adopt a more dynamical approach here. 

Equation (24) may be written in terms of the current density j, 

j 

v e = v , (25) 

nie 

where n» is the ion number density (taken equal to the electron's), and 
the final ion-neutral drift term is ignored. It is clear now why the Hall 
term on the right is important, with n{ small and the currents large 
enough to generate fields of interest. The induction equation (23) takes 
the form 

«_VX,. XB -^). (26) 

which introduces Coriolis-like terms in the behavior of the field lines. 
This makes the behavior of the fluid and the field symmetric: "field- 
freezing" tries to couple the magnetic field and fluid velocities, while 
Coriolis forces and Hall effect terms cause internal epicyclic motions in 
the matter and field respectively. Relative to the disk rotation, the fluid 
epicyclic motion is retrograde, while the sense of the field line motion 
depends on the magnetic helicity ft • B. 

A more clear picture of the effect of the Hall term emerges, if we 
consider its behavior in a uniformly rotating disk, ignoring for the 
moment any effects of finite electrical resistance. If we once again 
consider perturbations to a vertical field B = Be z , and denote linear 
perturbations by 5, plane wave disturbances of the form exp(ikz — iujf) 
satisfy the dynamical equations, 

ikB 
-wSvr - 2nSv (P - - — 5B R = 0, (27) 

ikB 
-iujdvs + 2VL5v R - - — SBs = 0, (28) 

47T/9 
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and the induction equations, 

k 2 Bc 
-icvSB R + 5Bs - ikB5v R = (29) 

-w5Bs - ^-^5B R - ikbbvA, = 0. (30) 

47renj 

Note the symmetrical appearance of 5B and 5v in the above, and that 
the Hall terms enter for the magnetic field precisely as the Coriolis 
terms do for the velocity field, except that they may induce a different 
sense of rotation. 

The dispersion relation for the above system is 



4 2 

id — CO 



2k 2 u 2 A + in 2 + kUH 



+ k A (u A + uf I ) 2 = (31) 



An 2 

where we have defined the Hall velocity uh by 

2 nBc , s 

"" - w (32) 

where c is the speed of light. Note that u 2 H may be positive or negative 
depending upon the sign of helicity factor Q • B. The Hall term alters 
both the effective epicyclic frequency of the coupled field-fluid system 
as well as the radial centrifugal force on a displaced fluid element, 
discussed in §2.2. The contribution of the magnetic field to the return 
force on a displaced element goes from k 2 u\ to k 2 u 2 A + k 2 u 2 H: and 
a simple measure of the importance of the Hall currents is the ratio 
x = u 2 H /u 2 A . 

The dispersion relation in the presence of differential rotation and 
resistivity is more complicated, and seems best left to the literature 
for discussion (e.g., Wardle 1999). But the stability requirements are 
easily stated. In the absence of resistivity one simply replaces u 2 A with 
u\ + u 2 H , 

oin 2 

k 2 u\(\ + x) + — — > STABILITY (33) 

am it 

In the presence of a resistivity r\ (units: cm sec -2 ), matters are a bit 
more complicated: 

dQ 2 
k 2 u 2 A [l + x + (/U + x)} + —— >0 STABILITY (34) 

d In R 

where 

C S ^. (35) 
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The effect of resistivity is to force the instability to longer wavelengths, 
lessening the perturbation currents. But once the wavelength exceeds 
the available linear dimensions, the system is stabilized. 

Overall, however, the effect of Hall currents on the stability of accre- 
tion disks are not simple. They both stabilize and destabilize, depend- 
ing upon whether ft is aligned or counter aligned with B. One relatively 
straightforward effect which may be seen either in equation (33) or (34) 
is that when B -» 0, i > 1 and the Hall term dominates. The same 
magnetic tension force is provided by much longer wavelengths than 
would have been possible without the Hall currents. This extends the 
range of unstable field strengths before dissipation processes kill the 
instability. 

We have seen that Hall effect currents affect the stability of a low ion- 
ization protostellar disks. How might they affect the nonlinear develop- 
ment of the turbulence? We are some ways from answering this, but one 
potential complication can be addressed. The Poynting contribution to 
the turbulent energy flux is 

Bx(v e xB). (36) 

A key piece of the alpha disk formulation is that the energy flux be 
a linear function of the same turbulent flow quantities responsible for 
angular momentum transport: the radial drift velocity and the R<p com- 
ponent of the stress tensor (Balbus & Papaloizou 1999). The Poynting 
contribution is one of the dominant components of the energy flux, and 
the fact that v e appears above instead of v looks like a complication — 
we have emphasized the importance of the difference v — v e . But here 
it is the azimuthal component of the unperturbed rotational velocity 
v e which is responsible for the radial contribution to the energy flux. 
(What is important for a perturbation may not be important in the 
unperturbed state.) For the Hall term to be negligible in this context, 
we find 

u\ < (m)(m) (37) 

where £ is defined by |V X B\ ~ B/£, a sort of Taylor microscale 
(e.g., Tennekes & Lumley 1972) for MHD turbulence. If £ is not much 
smaller than H 2 /R, this inequality will be amply satisfied, and the 
alpha formulation of MHD turbulence should remain intact. 
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4. The Role of Self-Gravity and Hydrodynamical Processes 

4.1. The Intermediate Disk 

Beyond 1 AU, the body of the disk is likely to be magnetically active. 
Rem drops rapidly until one reaches beyond 10 AU or so, where cosmic 
rays or X-rays can penetrate and leave a sufficient residual ionization 
to ensure good coupling once again. Although the focus of this meeting 
has been the zone of terrestrial planets, some mention of how the inner 
and outer disk are related is needed. 

It is now clear that many of the usual hydrodynamical mechanisms 
invoked in the past — turbulence generated by infall, convection, etc. — 
will not transport angular momentum outward. One that certainly does 
work, at least in principle, is self-gravity (e.g., Nelson et al. 1998). If 
the mass of the disk within a radius R is in excess of (H/R)Mq, it 
will become Jeans unstable, developing local (possibly global) spiral 
structure, which in turn is an effective outward transport mechanism. 
At 5 AU, the required interior disk mass is ~ 0.05 M@. This is much 
in excess of Jupiter's 1O _3 M0, but it is not unphysical. It is what 
one might guess if planet formation is efficient at the level of a few 
per cent, and consistent with the figure of 10 -2 M@ inferred for the 
accreted mass of the inner disk in FU Orionis outbursts (Hartmann, 
Kenyon, and Hartigan 1993). 

While this behavior may be common, gravitational replenishment 
is not likely to be universal, and certainly cannot go on indefinitely! 
What happens when the disk mass is too small? In some cases, we 
would expect to find passive disks with a hole on scales less than AU. 
At later stages in a disk's evolution, this will be unavoidable, as mass 
is lost one way or another. Strom, Edwards, and Skrutskie (1993) have 
found evidence for disks with inner holes on precisely these scales. They 
argue that these probably represent later stage transition disks when 
accretion begins to cease, that the inner parts of the disk need to 
be optically thin, and that the inner disk region is "isolated" from 
the outer. Masses of the optically thick inner disks were estimated at 
between 0.01 and 0.1 M & , consistent with self-gravity requirements. 
At the time of the Strom et al. review, planet formation was thought 
to be the most likely explanation. But among other possibilities (e.g., 
Konigl 1991), these disks may be candidates for magnetic depletion. An 
interesting prediction is that the size of the inner hole should correlate 
with disk metallicity in this picture: higher abundances should produce 
bigger holes, since the ionization levels needed for higher alphas could 
be maintained farther out. Detailed modeling is needed to test these 
striking general agreements of scale and predicted morphology. 
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4.2. Why Hydrodynamical Turbulence Fails 

The central feature of the picture of the Solar Nebula we have discussed 
here is that turbulent transport is regulated by magnetohydro dynam- 
ics, and that hydrodynamical transport is via self-gravity If differential 
rotation by itself were for some reason unstable, we would be forced 
to a completely different picture, one in which the turbulence prop- 
erties were much more uniform. This point of view retains prominent 
advocates (e.g., Dubrulle 1993, Richard & Zahn 1999), and should be 
addressed. 

The classical laboratory set-up for studying hydrodynamical turbu- 
lence in differential rotation is to examine flow between two rotating 
cylinders, the Taylor-Couette experiment. According to the classic text 
of Landau & Lifshitz (1959, p. 110), flow predicted to be stable by 
the Rayleigh criterion is shown by experiments to break down into 
fully developed turbulence at sufficiently high Reynolds numbers. The 
textbook by Zel'dovich, Rumaikin, and Sokoloff (1983, p. 321) is even 
more explicit on this point, arguing that the nonlinear breakdown of 
Keplerian flow is supported by experiments. 

These are extraordinary claims. We have found no such experiments 
in the literature; almost certainly none exist. The closest approxima- 
tion to a nonlinear breakdown of Keplerian flow is seen in experiments 
dominated by the rotation of the outer cylinder, or with a very small 
but positive angular momentum gradient. The first is characterized 
by (Kl/hidR 3> 2£l (shear dominating Coriolis forces), the second 
by a nearly vanishing value of n 2 (small vorticity). Though decid- 
edly nonlinear in their stability properties, each of these flows differs 
significantly from Keplerian, which has comparable values of k = f2, 
dQ/dlnR = —1.517, and 2Q. No Couette flow experiment has show any 
sign of instability when the inner and outer cylinders follow anything 
close to a Keplerian profile. 

To be fair, the Russian texts are somewhat misrepresentative; most 
proponents of hydrodynamical turbulence would concede that the ex- 
periments performed to date do not in fact show that Keplerian profiles 
are unstable. But, it is asserted, this is just a matter of going to higher 
Reynolds number, Re. Keplerian flow will be unstable at high Re, 
because that is the way shear flows behave. 

This intuition seems born of the classical bounded planar Couette 
and Poiseuille flow, and it presumes that the only dimensionless number 
characterizing the flow is Re. It also presumes an exquisite sensitivity to 
Re, since proponents are forced to argue that high resolution numerical 
simulations, which can easily recover the onset of laboratory Couette 
instabilities, once again simply have not achieved high enough Re. 
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To begin with, there are, of course, two nondimensional numbers 
characterizing Couette flow, both of which influence stability. One is 
the Re, the other is the Rossby number Ro, 



Ro 



ldlnQ 



2dlnR 



(38) 



which measures the importance of inertial forces relative to Coriolis 
forces. Obviously, Ro is a critical parameter in determining flow sta- 
bility. As noted explicitly in the review of Bayly, Orszag, & Herbert 
(1988), the effect of viscosity on unbounded flows (this includes shear 
layers and disks) is that of a regular perturbation, and the instabilities 
affecting such systems are essentially inviscid. The viscosity, when in- 
troduced, perturbs the flow at a level proportional to Re^ 1 <C 1, not, as 
in boundary layer theory, at the order unity level. Numerical Reynolds 
number effects are very small. The Rossby number is the only inviscid 
flow quantity of interest; to within a factor of 2 it is just the power law 
index of CI. 

The nonlinear breakdown of shear flow has been studied in some 
detail (Bayly et al. 1988), both quasi-analytically and fully numeri- 
cally. In all cases, it involves a two-stage process: an neutrally stable 
nonlinear two-dimensional equilibrium is first perturbed on top of the 
fundamental flow. This is a new "equilibrium." In the next stage, this 
new solution is itself linearly perturbed in three dimensions. The new 
equilibrium is exponentially unstable. The existence of neutral behavior 
in the intermeditate two-dimension solution is critical. It is a feature 
of systems whose linear restoring forces effectively vanish, as is the 
case for simple shear layers and constant angular momentum disks. 
It accounts for why nonlinear disturbances quickly become dominant. 
Systems similar to Keplerian disks, with comparable epicyclic and shear 
rates, do not allow the formation of neutrally stable nonlinear equilibria 
on top of the differential rotation. 

Finally, we note that turbulence maintains itself in shear layers by 
vortex stretching (Tennekes and Lumley 1972). Stretching is possible 
in neutrally stable flows for which initially neighboring elements drift 
apart when perturbed. This includes disks with nearly constant angular 
momentum profiles; these are indeed nonlinearly unstable. But in a Ke- 
plerian (or similar) disk, epicyclic motions keep elements localized, and 
vortex stretching is extremely inefficient (Hawley, Balbus, and Winters 
1999). 

Formal analytic arguments of the moment equations for fluctuations 
can be advanced to support these arguments (Balbus & Hawley 1998). 
The basic point is that in each of the unstable laboratory flows a 
dominant source term is present in their equation, which acts obviously 
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differently for flows in the stable regime. Numerical PPM simulations 
on high resolution grids (128 3 ) have been performed recently (Hawley 
et al. 1999) on local Keplerian flow. The effective Reynolds number 
at this resolution is very high, probably in excess of 10 5 . The analytic 
arguments were confirmed in detail, and not a trace of instability was 
observed. By way of contrast, our group has found that nonlinear insta- 
bilities in constant angular momentum flows can be recovered on grids 
as small as 8 3 ! When turbulence is directly imposed on a Keplerian 
disk in the form of driven convection, the angular momentum gradient 
compels inward angular momentum transport (Stone & Balbus 1996); 
when imposed by random dynamical forces, no transport at all occurs 
(Balbus & Hawley 1998). It is difficult to understand how supposedly 
diffusive numerics could be at once large enough to damp Keplerian in- 
stability on a very fine grid, small enough to allow nonlinear instabilities 
on an extremely crude grid, and small enough to allow imposed thermal 
fluctuations to drive inward angular momentum transport! Each of 
these behaviors is, however, a completely straightforward consequence 
of inviscid flow dynamics. It is simply untenable to cling to the notion 
that this is all somehow a gross misinterpretation of numerical artifacts. 
An angular momentum gradient in the flow is no less stabilizing than 
an entropy gradient in a stratified fluid, to which it is exactly analogous 
(Jeffreys 1928). One would not argue that the radiative interior of 
the sun is convectively unstable because the temperature gradient is 
a source of free energy that nonlinear disturbances can tap into. One 
would not argue that the Richardson criterion for stability, which is 
due to the presence of gravitational stratification across a shear layer, 
is something that disappears at large Reynolds number for nonlinear 
perturbations. Angular momentum stratification stabilizes shear flows 
as effectively as gravity stabilizes a bottom-heavy fluid. A century of 
investigation has produced no evidence to the contrary. 



5. An Overview of Solar Nebula Dynamics 

We are at an interesting and exciting stage in our developing under- 
standing of solar nebula dynamics. Nothing like a complete model is at 
hand of course, but even understanding how an isolated magnetized gas 
disk around a nascent sun behaves would represent important progress. 
And we have made progress. 

We have at hand an understanding of the physical basis for gen- 
erating and maintaining turbulence, and of enhancing transport in 
magnetized disks. The magnetorotational instability has broad analytic 
and numerical support. 
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We have not just an overall understanding of the way turbulence 
could work, we also have a handle on many of the details: ioniza- 
tion balance, Hall currents, nonlinear dynamics, for example. We know 
that metallicity is important, since the alkalis regulate magnetic cou- 
pling. There is, therefore, a direct connection between metallicity and 
gas dynamics. It is more usual, of course, for the metals to regulate 
thermodynamical properties. 

Models in which accretion, mediated by an external ionizing agent, 
proceeds through the the exposed layers of the disk have been proposed 
by Gammie (1996) and Glassgold et al. (1997). However, a detailed 
study of how the magnetorotational instability actually behaves for 
this geometry and resistivity has not yet been done. 

To summarize: the MHD picture of the (isolated) solar nebula is 
one in which turbulence is present within ~ 1 AU and beyond ~ 10 
AU. Gravitational instability may be present in the early stages on 
intermediate scales, maintaining a mass influx and replenishing the 
inner disk. When the disk mass falls to the point where this is no 
longer possible, the inner regions of the disk will be depleted, an effect 
that should manifest itself as a deficit at near IR wavelengths. This 
spectral feature has been observed, but an inner planet, or disruptive 
stellar magnetosphere (Konigl 1991) are also possible interpretations. 

It is a pleasure to thank C. Terquem for several useful comments, 
which lead to a significant improvement of this manuscript. This work 
has been supported by NASA grants NAG 5-7500 and NAG 5-3058. 
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